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Abstract 

We present the computation of the direct photon production cross-section in 
perturbative QCD to all orders in the limit of high partonic center-of-mass energy. 
We show how the high-energy resummation can be performed consistently in the 
presence of a collinear singularity in the final state, we compare our results to the 
fixed order NLO cross-section in MS scheme, and we provide predictions at NNLO 
and beyond. 



1 Introduction 

Perturbative QCD provides accurate theoretical predictions for hard processes at high- 
energy colliders. Logarithmic corrections to the lowest-order cross-sections can be system- 
atically computed in the region of large hard scale Q 2 , A 2 <C Q 2 ~ S, by a renormalization 
group approach which leads to the factorization theorem of mass singularities pQ. How- 
ever, the TeV energy range opens up the two scale region A 2 C Q 2 C S, where the usual 
perturbative expansion receives large contributions characterized by logarithms of the ra- 
tio x = Q 2 / S. In order to recover the accuracy of the perturbative results, logarithmically 
enhanced small- x contributions to the hard cross-sections, associated to multiple gluon 
emission, must be resummed to all orders. 

Prompt photon production [2] is a relevant processes for the study of hard interactions 
in high-energy collisions. For example, it is the most important reducible background for 



the H — > 77 signal in the light Higgs scenario [3]. A thorough understanding of this 
process in the small-x limit is thus relevant to make predictions for the LHC. 

Currently the direct photon cross-section is known up to 0(aa 2 s ) jl] and Sudakov 
resummation effects have been computed up to NLL accuracy [5]. Prompt photon pro- 
duction is especially useful to probe the gluon parton density over a wide range of x [6], 
since the initial state gluon appears already at leading order, and in particular the small- 
x (x > 1CT 5 ) region. For this reason we expect that, in prompt photon production, 
high-energy resummation affects significantly the fixed-order results. 

The general procedure for the small- a; leading-log (LLx) resummation of hard coef- 
ficient functions is well established in perturbative QCD within the framework of the 
fcj-factorization theorem [II IB], and involves the computation of the leading amplitude of 
the process with off-shell incoming gluons. This technique has been used to obtain re- 
summed cross-sections for heavy quarks photo- and hadro-production [3 [HJ EH ED] , deep 
inelastic scattering [Bl H2], Higgs production [TBI HU [T3] and recently for the Drell-Yan 
process [T6] . 

Prompt photon cross-section contains two different contributions: the direct compo- 
nent, where the photon participates at leading order to the hard process, and a frag- 
mentation component, which is needed to take account of the hadronic component of the 
photon. From a phenomenological point of view, at high-energy both terms are impor- 
tant [T7] ; in this work we will consider the direct contribution, leaving the fragmentation 
component to a future work. 

All the processes for which small-x resummation has been performed so far are free of 
collinear singularities in the final state since the corresponding cross-sections are totally 
inclusive; on the contrary such a divergence does appear in direct photon production 
because the process is exclusive with respect to the final state photon, which from this 
point of view must be viewed as another hadronic state [IB]. In this work we perform 
the high-energy resummation of the direct photon coefficient function consistently with 
the MS scheme of subtraction of the final state singularity to all orders in perturbation 
theory. 

2 Prompt photon production 

2.1 Collinear factorization 

The prompt photon process is characterized by a hard event involving the production of 
a single photon. Let us consider the hadronic process 
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According to perturbative QCD, the direct and the fragmentation component of the in- 
clusive cross-section at fixed transverse momentum q of the photon can be written as [5] 



q 3 ^ 7 ^ ±,q - = ^ / dxi f a/Hl {xi,fi 2 F ) f dx 2 fb/H 2 {x2^ 2 F ) x 

ab J X± Jx ± /X! 

x / dx<5 (x- — — J C2 b {x,a s {n 2 );q 2 ,n 2 F ,n 2 f ) + 

JO IV X\X2J 

+ ^/ dz z 2 d c/7 (z,jj 2 f )5 (x- -^T^j ^ b (x,a s (/i 2 );q 2 ,^,/iJ)| , (2) 

where we have introduced the customary scaling variable: 

4q 2 

x ± = ^-, 0<x ± <l. (3) 

in terms of the hadronic center-of-mass energy S = (P1+P2) 2 . In the factorization formula 
eq. (j5J) we have used the short-distance cross-sections 

rJ (c) 3 da ab ^ (c) (x, a a (// 2 ); q 2 , (j? f , h)) 

Lab ~ q dq ' 1 j 

where a, 6 and c are parton indices (g, q, g) while fi/H Hp) * s the P ar t° n density at the 
factorization scale fj,p. The fragmentation component is given in terms of a convolution 
with the fragmentation function d c / 1 (z, n 2 ). 
If we define the Mellin moments 



a f(N) = I dx ± x^~ 



1 ( 3 c/cr 7 (x_L,q' 



<iq 

Fi /H (N^ 2 F ) = / dx x N-1 xfi/ H {x^ 2 F ) 
Jo 



1 

_D c / 7 (iV, /i 2 ) = / dx x N_1 x 3 d c /ry(x, /i 2 ) 
Jo 

C^(iV) = fdxx N ^Ct\x) 
Jo 



the collinear factorization theorem in iV-space becomes 

a\N) = ^2F a/Hl (N)F b/H2 (N) (&JN) + J2 D ^( N )^( N ) ) ■ ( 5 ) 
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2.2 Leading-order coefficient functions and beyond 



At leading order the processes that contribute to the direct component of the prompt 
photon cross-section are 

qq^ig, q(q)g -> 7?(?)- ( 6 ) 

The corresponding leading-order hard coefficient functions are given by [21 HE] 



~1,LO/\ q 3 ^g(g)g->7g(g) ^ 2 * 



X ^ X s 



in terms of the partonic variable x = 4q 2 / s, where s is the partonic center-of-mass energy. 
QCD corrections to the Born coefficient functions eqs. (CD) and (jSJ) have been computed 
in ref. jl] up NLO accuracy. 

In the high-energy limit, the leading coefficient function is logarithmically enhanced 
by contributions of the form 

oo 

W = C l L °W + ««2 E c S } («slogx) fc + O {aa 3 s (a s \ogxr) (9) 

fc=0 

which, in A^Mellin space, becomes a sum of poles at N = of increasing order: 

OO £ 

q 9 (iv) = q^) + ^E^ } (|) +°H(!) n )' n — (io) 

As in the case of heavy quark production (HQ) and of the Drell-Yan processes (DY), the 
Born coefficients C 1 ^ are regular as N — > 0: indeed, eqs. ( ITfBi) vanish when x approaches 
zero. The first singular term in eq. ( flOj) is a simple pole in N = given by the NLO 
contribution to the perturbative series; in x-space, this pole corresponds to a constant 
value, while the small-x logarithms arise from the poles of increasingly higher order. 

The NLO contribution to the expansion given in eq. © aa^cfj has been computed for 
various processes in ref. [IH] , in particular for direct photon production, by considering the 
Feynman diagrams in fig. (Op) where an extra gluon is radiated from the initial state. This 
is a general feature of the high-energy limit to all orders in perturbation thory: dominant 
contributions at high-energy (LLx) are given by the exchange of spin-1 particles in the 
t-channel therefore all the relevant Feynman diagrams in the small-x limit are given by 
the BFKL ladders in fig. (El). 

From the NLO onwards, the direct photon cross-section acquires a final state diver- 
gence when the photon becomes collinear to the outgoing parton; this collinear divergence 
cannot be removed by adding the virtual corrections, rather it must be absorbed in the 
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Figure 1: In this picture are shown the dominant Feynman graphs in the high-energy 
limit. Both of the two diagram determine the constant behaviour of the hard coefficient 
function Cgg at small-x. 

fragmentation component of eq. (j2J) as it happens for the initial state collinear divergences 
which are properly absorbed in the definition of the parton densities [18]. In the next 
sections we will show how to remove this divergence in the MS scheme consistently with 
the resummation procedure of the gluon ladder and we will give the analytic expression 
for the resummed hard coefficient function Cgg and Cgg. 

3 High-energy resummation 

In this section we compute the leading logarithmic corrections to the direct photon cross- 
section at high-energy. The high-energy resummation of gluon ladders arises from the 
general formalism of the fct-factorization theorem [7J [8] and can be performed by following 
the computational procedure outlined in ref. [Hj which allows us to compute the coefficient 
of the dominant log x to all orders in perturbation theory. This procedure requires the 
calculation of the Born cross-section with off-shell (transverse) incoming gluons and off- 
shellness fixed in terms of their transverse momenta kj (impact factor). 

For a single off-shell gluon, we can parametrize the dependence on the virtuality k 2 
through the dimensionless variable 

e = £. (id 

where Q 2 is the hard scale of the process which determines the argument of the run- 
ning coupling a(Q 2 ). In direct photon production Q 2 is the magnitude of the transverse 
momentum of the photon q 2 . The LLx resummation is performed by taking a Mellin 
transform of the off-shell cross-section 

h{N,M)= da M ~ l dx x N ~ 1 a{x,i), (12) 
Jo Jo 
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Figure 2: Multiple gluon emission (BFKL ladder) from the initial state in direct photon 
production. 




Figure 3: Diagrams for direct photon production. 



and by identifying M as the sum of leading singularities of the largest eigenvalue of the 
singlet anomalous dimension matrix [20] (BFKL anomalous dimension) 

M - m) + °© < 13 > 

"?•(!) = f>(^rf. c„ = l,0,0,2<(3). (14) 

n=l ^ ' 

This corresponds to the sum of the high energy contributions coming from all diagrams of 
fig. (j2). Finally, the N fc LO coefficient of the maximum power of logx is given by expanding 
the impact factor in powers of a s . 
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3.1 The off-shell cross-section 



Let us consider the process 

9*{k) + q{p)^j{q) + q{p'), (15) 

with an off-shell incoming gluon. We use a Sudakov parametrization for both incoming 
and outgoing momenta of the diagrams in fig. [31 thus we have 

V = Z2P2, (16) 
k = zipi + k ± , (17) 
q = Zi^iPi + (1 - x 2 )z 2 p 2 + q±, (18) 

V = (1 - Xi)ziPi + X2Z2P2 + k± - q±, (19) 

where 

^J_ ~~ I I ~~ k 1 

ql = ~\q±\ 2 = -q 2 , 

with light-like vectors p\ and p 2 such that p\ ■ p 2 = S/2 where S is the energy in the 
center-of-mass frame. The relevant scalar products are 

p ■ k = s/2, 
p ■ q = X\s/2, 

q 2 - 

k-q = (1 -x 2 )s2 - k± ■ qj_ = k ± -q ± , 

2x\ 

where we introduced the longitudinal energy s = Z\Z 2 S. The ci-dimensional phase space 

d$(<*) = - dVV 5G» /2 )5(g 2 )5 (d) (k + p-q-p') = 

(27r) d 2 

rfV(b+fc-?) 2 )% 2 ) (20) 



(2tt 



id— 2 



can be rewritten in terms of the Sudakov parameters since d d q = ^d d ~ 2 c\dxidx 2 . We 
obtain 

1 s 

= ^- 2 q rf XlC / X2 5(_k 2 + s-XiS- (1 -x 2 )s + 2k-q) • 

(27r) 2 2 

■<5(-q 2 + Xi(l - x 2 )s) = 

- yg) ^"-^*-^-?^^ (2i » 

where # is the angle between q and k and we used the last 5-function to perform the 

2 

integration in x 2 , now fixed to x 2 = 1 — jL, which implies x < 4xi. 
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Since we are interested in the differential cross-section a(x,£) = q 3 ^ 2 -, our phase space 



in four dimension reduces to 
,d$ (4) . 1 



q 



dq 



(2vr) 2 2x 



q dv o 



s(l - si) - 



2k-q 0(s-k 2 ), (22) 



which in terms of the dimensionless partonic variables x = 4q 2 /s and £ = k 2 /q 2 reads 



1 XSC?Xi 



(2tt) 2 
1 



d8 5 



iXi 

sxdxi 



4 



.r 



(1-zi) 



C - - + 2v^cos^ 0(- - 0e(4x x 

Xi J X 



X) 



(2tt) 2 32xi(1 -xi] 



1 £ + ^-2^cos# 



.r 



4(1 -xi) 




0©(4xi -x). 



In the resummation procedure of refs. El E] the computation of Feynman diagrams 
is performed by using the eikonal rule for the gluon polarization sum 



~k 2_; 



(0,k,0), 



(23) 



understood as the projector V over the high-energy singularities, analogously to the ap- 
proach of refs. [21], which factorizes the gluon ladder from the Born coefficient. The 
channels s and t lead to the simple result for the amplitude in d = 4 + 2e dimensions 

— — Ap 2 n 2 V( ~ 2 

a^(x, Xi ,o = y.m 2 - Ae9s 



4e 2 # 



2-2C A 



(q z 



sxi) 2 + s 2 x\ + eq 4 
sx\ (s 



k 2 ^ 



2^2 



lQe 2 g 
2-2C A 



1 - ~T~ 



^i(s-0 



(24) 



averaged over color and helicity (of the quark) and summed over the final states. As 
shown in eq. ffTUl) . the high-energy enhancement appears as a series of poles in N = 0, 
therefore we are interested in the most singular term in the small- N limit. 

Since the off-shell cross-section is well behaved at N = 0, all the singular terms come 
from the substitution shown in eq. (TIBl) . hence at this level we can reduce the computation 
to the N = moment of the impact factor in the (N, M) space 



h(0,M) 



1 



(2tt) 2 2s 

f4/x 



2 77 



dO 



s dxi 



32xi(l — xi 



Li: i 



dx ■ 



o 



dU M ~ 1 A^(x,x 1 ,0S - 

x 



^ + ±-2^cos6 



(2 



1 1 
~2 



7T 



2k 



+oo 



d0 

dp 



dxi 



4(1- 

M-l 



32xi(l — xi j 
2 AV>(l/p, Xl ,Z)6 



d^ 



e+i-2v^cos^ 



4(1 - Xl ) 



(25) 



S 



where we introduced the variable p — 1/x and we exchanged the order of integration of £ 
and p. By using the delta function we obtain: 



(27r) 2 2 7 ' J 32x 1 (l-x 1 ) 



(26) 



where we have defined 

P= 4(1 -x.) ' (27) 
The argument of the Heaviside O-function in eq. (1261) is always positive since 

1 _ ex 1 + l-2x lv ^cosg-l + x 1 _ £-2v^cosg + l , (yg- l) 2 
^ 4xi 4xi(l — xi) 4xi(l — Xi) 4xi(l— Xi) 

when 1/xi > £, and 

£ £ + ^_2 v ^cos^-£(l-x 1 ) ^-2 v ^cos^ + £x 1 (-^-v^) 2 n 

p— - = — = — > > 

P 4 4(1 -xt) 4(1 -xO 4(1 -a*) 

in the opposite case. Therefore we have 

/i(0, M) = ^I^ dff -^^^ I ^_ y J J ^ (1M . 1 , f) . (28) 



The integration over the region < 9 < 2n, < x\ < 1 in eq. ( |28|) is always divergent 
when £ > 1, z.e. |k| > |q|; indeed the latter condition defines the kinematical region where 
the photon can be radiated collinearly to the quark in the final state. In the collinear 
limit, the amplitude in eq. (1241) is singular and the divergence is given by the fermionic 
propagator in the s-channel. In the collinear limit we have: 

£x 2 - 2^0030 + 1 ^ 
Xl — xf 



which happens when 

# = 

Xl = 7f 



(30) 



9 



3.2 Subtraction of the collinear singularity in MS scheme 



In order to cancel the collinear divergence in the MS scheme, first, we regularize the 
integrations of eq. ( !28l) by subtracting the collinear limit of the four dimensional amplitude 
before doing any integration, second we recover the pole in e = and the remaining finite 
parts by adding back the same quantity computed in d dimensions. 
We can do this by writing the impact factor in 4 + 2e dimensions as 

= / dct) {d) A {d) . (31) 



We then remove the singularity of the integrand by introducing a function T>^ which has 
the same singular behaviour of the squared amplitute A^ in the collinear limit eq. flHUj) . 
In four dimension we have 

^ = ^(i/VI) (1 _^ + , 2 eK-i), (32) 

where 

P[z) = I+U^£)! (33) 

z 

By adding and subtracting the phase space integral of the function T>^ to the <i-dimensional 
impact factor we obtain 

h {d) = lim (J d(f) id) A w - J d^V^A + j d^V^ + 0(e) = 

= J rf0 (4) {A {A) - P (4) ) + Ia + J d(f) {d) V^ + 0(e) (34) 

where the first integral is finite in four dimensions and the finite part fa comes from the 
linear term in e in the <i-dimensional amplitude eq. (j24p . 

By using the <i-dimensional phase space, the last term in eq. ( J341 is 



~2C A T(l/2 + e) U 2 J £ ^ NQ [{I ~ Vlxi) 2 + 2 ] ' { ) 

where the dimensional scale fi 2 (introduced by dimensional regularization) from now on 
will be identified with q 2 . By using this result in eq. (I34p and performing the Mellin 
integrations with JV = 0we have 

aa , ^(4.)-. i r dxi p j^nfMi-D 



2C A T(l/2 + t)iJ l J_e [(l-v^i,) 2 + » 2 ] 
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where we have extended the limits of the angular integration of the function T> in order 
to simplify the results of the integration while the finite part f\ is 



naas 1 ( 7 , 

~ C A 2-Af [ 

All the integrations in eq. (13"6"1) can be performed in closed form, thus by subtracting 
the pole in e = with the usual combination 1/e + 7^ — log47r we obtain^ 

hm = ^{e( £ -i)(-i + ^(i-u*2)) + 

+ sign« - 1) (3 - (l + 1) log(l - ?) 2 - ^ log (^) 2 ) } , (38) 

while the M-Mellin moments are 

» ^ ttct s i f (7 - 7M + 2M 2 ) / . . TT 2 \ 1 

M °' M) = OT { (M-l)(M-2)(2M-3) ( TC ° t(Ml) + 2H "-> + M^) + 2^M 

(39) 

where Hm-2 is the harmonic number of argument M — 2. 

4 Results 



4.1 The resummed coefficient function 



By expanding the impact factor obtained in the previous section around M = 0, we obtain 
ww w x 7raa s f 7 67, r 385M 2 /2323 7C(3)\ 

/ 14233 49am /87307 331«3) 7£(5)^ M * + (M «) \ . (40) 

V 7776 18 J V 46656 108 3 / J 



^dn ref. [19j the impact factor h q (a) was only computed in the region < £ < 1 where no collinear 
singularity is present. The result for h q , given in eq. (3.7) of ref. |19j . is seen to agree with our result 
recalling that in the region < £ < 1 h q is related to h eq. (|3"5|) by 

f^P = M0, (o<e<D. 

Note however that the expression for the cross-section er,j ff (q > pt) given in eq. (3.11) of ref. [115] is too 
large by a factor 2 (22] . 
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Notice that in this formalism the collinear divergence from the initial state appears as a 
simple pole in M = 0. The resummed coefficient function in the MS factorization scheme 
is given by the relation 

q g (N,a s ) = Mh(N,M;a s )R(M)\ M=ys (41) 
in terms of the impact factor eq. ( l39l) and the function 



R(M) = 1 + ^CsM 3 - jjc 4 M 4 + 0(M 5 ) (42) 

which takes into account finite parts coming from the MS subtraction of initial state 
collinear singularities and where M is identified as the BFKL anomalous dimension 

7,m=^ + 2C s m 4 + 2 C! mY.., = (43) 

We have 



NJ N ^ \NJ ™ \NJ 7T 



~ 7 iraa s f 7 67 (a s \ 385 ( a s ^ /2323 49C(3)\ / 



a N 3 



a,\ 4 „ //a N 5 



/ 14233 _!^ + 308C(3)' 

V 7776 720 27 y Viv/ \\iv/ y / 1 ; 

The NLO term in eq. (|44p gives, in the x-space, the constant value 67/36 ao? s which is in 
agreement with the fixed order calculation of refs. [HE]. By using the high-energy color 
charge relation between the hard coefficient functions 



A 

we can also obtain the LLx contributions coming from the process qq — > 



~ 7 qgC F [67 385 ( a s \ ^2323 49((3) ^ /e* s N 2 

««) = « ]v G 1 36 + 216 liV ) + (,1296 + J \N ) + 

/14233 _ 77T 4 308C(3) \ /o^ 3 //o«\ 



V 7776 720 ' 27 J ViV ) + ° ^ 



4.2 Phenomenology 

In fig.|4]we compare in x-space the coefficient function C qg (x, a s ) at LO, NLO and N 4 LO in 
the high-energy limit. The large contributions at small- x spoil the perturbative expansion 
and must be resummed in order to recover accurate results. The resummation of these 



12 




Figure 4: The picture shows a comparison between the leading order direct photon coeffi- 
cient function (black solid) and the relative small- x corrections up to NLO (green dashed) 
and N 4 LO (red dashed). 



logarithmic terms in the hadronic cross-section can be performed to all orders [121 [23] 
including running coupling. A full phenomenological study could be performed by com- 
bining the resummed hard cross section computed here with the resummation of GLAP 
evolution equations, following the formalism of refs [121 [23] (see also ref. [21] for an al- 
ternative approach). However, we can get a feeling for the size of resummation effects by 
comparing the result in eq. (|44p with the DIS coefficient function Cf(iV). As shown in 
fig. Othe ratio between Cf and C^ g (both of them normalized to the respective LO values) 
is of order 1, therefore we expect that, at low-x, resummation effects may be as important 
as those obtained in the DIS case. 



5 Summary 

In this work we performed the computation of the hard coefficient function for direct 
photon production in the high center-of-mass energy limit (the small-re limit) to all or- 
ders in perturbative QCD. We have shown how to resum the leading logarithms of the 
Bjorken variable x which stem from the BFKL gluon ladders, and how the high-energy 
resummation technique can be joined with the MS scheme of subtraction of a final state 
collinear singularities. 

High-energy resummation for both PDFs evolution and coefficient functions are well 
established tools which are needed to quantify the small- £ effects on perturbative results. 
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X 



Figure 5: Ratio between F 2 (Cms an d direct photon (DPhP) (C 7 ) coefficient functions, 
both of them normalized to their LO values 

A clear understanding of this high-energy effects will be crucial expecially for the LHC 
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